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O . Abstract 

We have performed stacking image analyses of galaxies over the Galactic extinction map constructed 
by Schlegel, Finkbeiner & Davis (1998). We select ~ 10 7 galaxies in total from the Sloan Digital Sky 
Survey (SDSS) DR7 photometric catalog. We detect clear signatures of the enhancement of the extinction 
in r-band, AA r , around galaxies, indicating that the extinction map is contaminated by their FIR (far 
, infrared) emission. The average amplitude of the contamination per galaxy is well fitted to AA r (m r ) = 

0.64 x io 017 ' 18_ " ir ' [mmag]. While this value is very small, it is directly associated with galaxies and may 
have a systematic effect on galaxy statistics. Indeed this correlated contamination leads to a relatively 
large anomaly of galaxy surface number densities against the SFD extinction Asfb discovered by Yahata 
Q-f et al. (2007). We model the radial profiles of stacked galaxy images, and find that the FIR signal around 

q . each galaxy does not originate from the central galaxy alone, but is dominated by the contributions of 

nearby galaxies via galaxy angular clustering. The separation of the single galaxy and the clustering terms 
c/3 . enables us to infer the statistical relation of the FIR and r-band fluxes of galaxies and also to probe the 

flux-weighted cross-correlation of galaxies, down to the magnitudes that are difficult to probe directly for 
individual objects. We repeat the same stacking analysis for SDSS DR6 photometric quasars and discovered 
— , ■ the similar signatures but with weaker amplitudes. The implications of the present results for galaxy and 

^ ' quasar statistics and for correction to the Galactic extinction map are briefly discussed. 

C . Key words: ISM: dust, extinction — cosmology: observations 
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| 1. Introduction the value of the extinction ^4 rj sFD for low extinction re- 

■ gions (A rj sFD < 0.1); see their Figure 4. After several 

t-H ' Galactic extinction is the most fundamental correction careful analyses, they concluded that this anomaly origi- 

. that should be applied to virtually all astronomical data, nates from the far infrared (FIR) emission from the SDSS 

The most widely used map for the purpose is constructed galaxies themselves. If such additional FIR emissions from 

by Schlegel, Finkbeiner & Davis (1998; hereafter, SFD). extra-galactic objects are incorrectly interpreted as those 

They first made dust temperature and emissivity maps from the Galactic dust, the extinction of those regions is 

K> ! (0.° 7 FWHM spatial resolution) from COBE/DIRBE data overestimated. The amount of the FIR emission would 

at 100 /im and 240 fjjn. Then a finer resolution map was be an increasing function of the surface number density 

created for dust emission (6.'1 FWHM angular resolution) of galaxies, and the correlation between them should be 

from IRAS/ISSA data at 100 using the COBE tern- visible especially where the Galactic extinction is small, 

perature map as a calibrator. Finally they constructed Moreover the correlation should become even stronger af- 

the maps of reddening and extinction assuming a simple tcr correction using the contaminated extinction value, 
linear relation between far infrared (FIR) flux at 100/im Although the proposal of Y07 for the origin of the 

and dust column density with a temperature correction for anomaly is reasonable, it is admittedly based on indi- 

dust emissivity. As is clear from their construction pro- rect and circumstantial evidence. In the present paper, 

cedure, the SFD map does not correspond to absorption- we report on a detection of FIR emission signatures from 

weighted, as required for most purposed in astronomy, but SDSS galaxies and quasars through a stacking analysis 

is an emission-weighted dust extinction map. of the SFD map. The amount of the FIR emission from 

This motivated Yahata et al.(2007; hereafter Y07) to those objects is statistically consistent with that required 

examine the validity of the assumed proportionality be- to explain the anomaly discovered by Y07. 
tween extinction and emissivity using the surface density The present paper is organized as follows; in §2 we de- 

of SDSS (York et al. 2000) galaxies in a given appar- scribe part of the SDSS DR7 data catalog that we use 

ent magnitude range as a calibrator. Surprisingly, Y07 in the analysis. After presenting the stacking analysis 

found that the galaxy surface densities increase against method, we show the stacked images and the resulting 
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fit to the analytical profile. Section 3 discusses the phys- 
ical interpretation of the extracted parameters for FIR 
emission from SDSS galaxies and comparison with the 
cross-correlation with IRAS PSCz catalog. We also re- 
peat the similar stacking analysis for SDSS DR6 photo- 
metric quasars in Section 4. Section 5 summarizes the 
conclusions. Finally in Appendix we show the IRAS Point 
Spread Function that we measure. 

2. Stacking analysis 

2.1. SDSS DR7 Data 

In the following analysis, we use the SDSS DR7 pho- 
tometric galaxy catalog, which covers 11663 deg 2 of sky 
area, with photometry in five passbands; it, g, r, i, and z 
(For more details of the photometric data, see Stoughton 
et al. 2002; Gunn ct al. 1998; Gunn ct al. 2006; Fukugita 
et al. 1996; Hogg et al. 2001; Ivezic et al. 2004; Smith et al. 
2002; Tucker et al. 2006; Padmanabhan et al. 2008; Pier et 
al. 2003). According to photometry processing flags, we 
carefully remove bad photometry data and fast-moving 
objects, which are likely in the Solar system or associated 
with the interplanetary dust. We also exclude objects in 
masked regions to avoid unreliable photometry data. For 
more details of our data selection, see Yahata et al. (2007). 

Our analysis below does not exclude the galaxies that 
are also detected by IRAS PSCz (Saunders et al. 2000; see 
§3) that are removed in the SFD map. We made sure that 
this has a negligibly small effect on our result because of 
the small number of overlapped galaxies. 

2.2. Stacking Method 

As discussed by Y07, the amount of FIR emission from 
SDSS galaxies that explains the anomaly is very small, 
and it is impossible to detect for individual galaxies. 
Therefore we stack those regions of the SFD map cen- 
tered at the positions of SDSS photometric galaxies (in 
this section) or quasars (in appendix) over their appropri- 
ate magnitude bins. For that purpose, we use a contiguous 
region in the SDSS DR7 photometric catalog as shown in 
Figure 1. 

The original SFD map divides all sky area into 2.'37 x 
2.'37 pixels and the extinction value is provided for the 
central position of each pixel. The histograms of A ri spD 
evaluated at those pixels as a function of the number of 
galaxies with 15.5 < m r < 20.5 within the pixel, -/Vg jP i x , are 
shown in Figure 2. While the overall shapes of the his- 
tograms are very similar for different Ag !P j x , the normal- 
ized probability density function (PDF) plotted in Figure 
3 exhibits the small but systematic shift toward the larger 
^4r,SFD with increasing iVg iP i X . This indicates the correla- 
tion of the Galactic extinction and the background galax- 
ies that will be extensively discussed in the present paper. 

First we show the result of stacked SFD map images 
centered at photometric galaxies in the r-band magnitude 
range of 17.5 < m r < 19.4 randomly selected from a con- 
tiguous region in Figure 1. In this procedure we evaluate 
the value of A-.sfd on 0.'2 x 0.'2 pixels over 40' x 40' im- 
ages by cloud-in-cell interpolation of the 4 nearest neigh- 
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Fig. 1. The region of the sky used for the present analysis, 
which is shaded according to the extinction value A r sFD- The 
yellow lines indicate the inner regions used for comparison in 
Subsection 2.3. 
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Fig. 2. Histograms of A r sfd f° r 2f37 X 2.'37 pixels over our 
selected survey region of SDSS DR7 as a function of the num- 
ber of galaxies within the pixel, iV g p ; x ; Af gjP ; x = — 3 in black, 
4 — 7 in red, 8 — 11 in blue, and 12 — 15 in green. 



bors in the original SFD pixels. Upper panels of Figure 4 
clearly show the presence of the strong feature of A rj sFD 
around SDSS galaxies, which becomes more pronounced 
as increasing the number of stacked galaxies. For refer- 
ence, lower panels of Figure 4 show the stacked SFD map 
images centered at the same number of randomly chosen 
positions from the same region of the corresponding top 
panels. This result directly confirms the interpretation 
of Y07 that the SFD map is contaminated by the FIR 
emission from SDSS galaxies. 

2.3. Radial profiles of galaxies 

To estimate the dependence of the contribution to 
^r.SFD, or cquivalently the amount of the FIR emission, 
on galaxy r-band magnitudes, we stack the images at the 



Detection of FIR Emission from SDSS galaxies 



3 




: ■ 79.4478 > 




■ - ■ 79.2779 



N,andom=1 00000 



4% 



Fig. 4. Stacked images of the SFD map for 40' X 40' regions; Upper panels are centered at the positions of SDSS galaxies of 
17.5 < m r < 19.4, and lower panels show the reference images centered at randomly selected positions. Left, middle and right panels 
correspond to images stacking 10 3 ,10 4 ,10 5 images. 
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Fig. 3. The normalized probability density function (PDF) 
of A r gpu corresponding to Figure 2. 

location of galaxies according to their r-band magnitudes. 
The results are plotted in Figure 5. Thanks to the sig- 
nificantly large number of the SDSS galaxies, those im- 
ages are highly circular, assuring that the signals do not 
originate from the Galactic contamination. All the angu- 
lar radii of the images are very similar to the expected 
smoothing length of the SFD map (= 6'.1 FWHM), and 
one may naively interpret that the central signal is domi- 
nated by the single galaxy contribution. As we will show 
below, however, this is not the case; in most cases the 
signal is rather dominated by the contribution from the 
nearby galaxies. Qualitatively this is understood from 
Figure 3 that shows the systematic increase of A ri sFD as 



a function of numbers of galaxies in the pixel whose size is 
much smaller than the overall smoothing size of the SFD 
map. 

To proceed more quantitatively, we attempt to model 
the radial profile of the stacked images as follows. Denote 
the angular profile of a single galaxy with r-band mag- 
nitude m r as E|(0;m r ). Then the statistically averaged 
profile of the stacked images centered at the galaxy is given 
by 



£* ot (0; mr) = K(0; mr) + K(9; m r ) + C, 



(1) 



where £g(0;m r ) denotes the clustering term correspond- 
ing to the contribution from the nearby galaxies, and C 
represents the background level of the extinction. Naively, 
C is expected to be independent of m r and computed from 
the PDF of the extinction P(A) (see Figure 3) as 



C=(A)± 



OA 



(2) 



where N g is the number of stacked galaxy images, and the 
mean and rms are given by 



(A)= / AP(A)dA, 
Jo 

/>oo 

a\ = \ A 2 P(A)dA- 
Jo 



(Ay 



(3) 
(4) 



As we see below, however, this is not the case. Therefore 
we treat C as a free parameter for each magnitude bin in 
the fitting analysis described below. Figure 6 plots (A) 
with quoted error-bars of a a as a function of N giPlx . 
The clustering term is written as 



E^(6/;m r ) = JJ dm'dip E|(0 - <p;m') 

dJV g (m') 



x w g (ip;m ,m r )- 



dm' 



(5) 
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Fig. 5. Stacked images of the SFD map for 40' X 40' centered at SDSS galaxies of different r-band magnitudes (m r = 15.5 ~ 20.5mag.) 
in 0.5 magnitude bin. The magnitude range and the number of galaxies in the range are denoted in each panel. 



where w g (cp; to', m r ) is the angular galaxy cross- 
correlation function between magnitudes to' and m r , and 
dN g (m')/dm' is the differential galaxy number density. 

Given the large smoothing length of the SFD map (6.'f 
FWHM), a single galaxy profile is expected to be approx- 
imated by the circular Point Spread Function (PSF), in- 
dependently of its intrinsic profile. Thus we adopt the 
Gaussian PSF profile: 



(6) 



£^(0;m r ) = I^(TO r )exp(-— ), 



where a is the Gaussian width of the PSF. The Gaussian 
approximation of the PSF is justified in Appendix. Also 
we assume that the angular cross-correlation function is 
given as 

w g ((p;m' ,m r ) = if (m',m r )(<p/<po)~ 7 , (7) 
where the constants ipo and 7 are assumed to be inde- 



pendent of to' and m r . We adopt 7 = 0.75 (Connolly et 
al. 2002; Scranton ct al. 2002), which is valid for ip < 1°. 
With equations (6) and (7) , equation (5) reduces to 



-^2 



x 1.F1 1 



;1; 



2a 2 r 



(8) 



where ii 7 i(a;/3;x) is the confluent hypergeometric func- 
tion, and 



E= u (m r ) =2na 2 



V2a 



r 1- 



7 



x / dm'Xf (m')K(m',m r ) 



dN s (m') 
dm' 



(9) 



Equation (8) results in the extended tail due to the clus- 
tering term in addition to the Gaussian tail of the single 
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Fig. 6. Mean values of A ri sFD of pixels containing iV g p i x 
plotted against A^ gjP i x . The quoted error bars indicate the 
corresponding rms. 

central galaxy. The latter is negligible at ^> a, and the 
observed tail of the profile around galaxies is basically 
dominated by the clustering term. 

The average radial profiles of the stacked images cen- 
tered at photometric galaxies are plotted in Figure 7. 
Filled circles and triangles correspond to galaxies in the 
different r-band magnitude ranges in Figure 5, and the 
quoted error-bars represent rms in each circular bin of 
A6 = 0'33. The solid curves are the best-fit model of equa- 
tions (I), (6), and (8). In fitting those curves, we treat 
Eg , E g °, and C for each magnitude bin, and the width of 
PSF, a as free fitting parameters. Here we assume that 
a is independent of m r , and fit it simultaneously for all 
profiles with different magnitudes. The resulting best-fit 
value of a = 3'1 is reasonable, given the resolution of the 
SFD map (2f 59 in Gaussian width) and the additional 
smoothing due to the 2. '37 pixelisation and our cloud-in- 
cell interpolation. 

The best-fit parameters for S g °(m r ) and £ g °(m r ) are 
plotted in Figure 8. Even at the central position of the 
stacked images, the FIR signals are indeed dominated by 
the clustering term E g rather than the single galaxy term 
(see Figure 9). Since the value of 7 is somewhat uncer- 
tain and also depends on magnitudes, m r and m' r , we 
performed the same fitting analysis by varying the value 
for 0.65 < 7 < 0.85, which roughly covers the range of 7 for 
our sample magnitudes. The dashed lines in Figure 8 show 
best-fit values for 7 = 0.65, and dotted lines for 7 = 0.85. 
Although the best-fit values of E g °(m r ) are affected by 
the choice of 7, especially for the faint magnitudes, the 
results for E g °(m r ) hardly change and still dominate the 
single term. 

In reality, we did not anticipate this before performing 
the radial fits. Therefore we directly measure the angular 
cross-correlation of SDSS galaxies for comparison, and in- 
dependently estimate the clustering term using equation 
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Fig. 7. Radial profiles of stacked galaxy images correspond- 
ing to Figure 5. Solid curves indicate the best-fit model of 
equation (1), (6), and (8). 

(5), which we call £ g ° (model). Dashed and Solid curves in 
Figure 10 indicate E g °(modcl) and £ g °(fit)+E g (model), 
respectively. We still adopt E g °(fit) (blue crosses), ob- 
tained by the profile-fits above, for the single term. The 
results indicate that resolved SDSS galaxies explain only 
25~40% of the clustering term. This implies that the re- 
lation between r-band and FIR luminosity has significant 
scatter and/or substantial fraction of FIR sources are not 
resolved by the SDSS optical selection criteria. The es- 
timated clustering term exceeds the single term by more 
than a factor of a few for faint magnitudes. 

The fitted values of the background noise term C are 
plotted against m r in Figure 11. According to our model, 
they should not depend on m r , but it is not the case at all. 
A systematic decrease of C against m r is clearly seen. We 
are not yet able to understand this systematic behavior. 
We repeated the same analysis by selecting those galaxies 
located in the inner contiguous regions (160° < a < 220° , 
5° < 5 < 80°; see Figure 1). The results are plotted in red 
crosses after shifting 25 mmag, just for the ease of visual 
comparison. The fact that the value of C is also sensitive 
to the region of the map and the weaker m r -dependence 
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Fig. 8. Best-fit parameters characterizing the FIR emission 
of galaxies against their r-band magnitude. Crosses, triangles, 
and circles indicate the best-fit value of S|° (fit), £g° (fit), and 
E g °(fit) +E g °(fit), respectively, assuming 7 = 0.75. Shaded 
regions indicate the variance of the best-fit parameters for 
0.65 < 7 < 0.85. Dashed curve indicates the results for 7 = 0.65 
and dotted curve is for 7 = 0.85. 
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Fig. 9. Ratio of the clustering term and the central galaxy 
contribution as a function of the r-band magnitude of the 
central galaxy. 

in this case may imply that this is due to a subtle selec- 
tion criterion difference of SDSS photometric galaxies as 
a function of the Galactic extinction itself, including the 
star-galaxy separation in different area, for instance. We 
do not investigate this issue in the present paper since 
our main finding is unlikely affected by the position de- 
pendent value of the background noise term. Actually, the 
best-fit values for other quantities are not sensitive to the 
particular choice of subregions in the SFD map. 

Incidentally the small value of C with respect to the 



Fig. 10. The clustering term Eg (model) estimated from the 
SDSS angular cross-correlation is plotted as a function of the 
r-band magnitude of the central galaxy (dashed curve). Solid 
curve indicates the sum of the single and clustering term 
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Fig. 11. The background noise level C against the r-band 
magnitude of the central galaxy. Black crosses indicate the 
results for all our sample and red ones are restricted in the 
inner regions shown as yellow lines in Figure 1. 

general trend at 16.0 < m r < 16.5 is the reason why the 
corresponding profile in Figure 7 does not follow the sys- 
tematic trend of the other profiles. 

We note that, if a is treated as a free parameter sepa- 
rately for different magnitudes, <r(m r ), the resulting fits 
in the cases of m r > 19.0 return a small negative value 
for We found that this is simply due to the fact the 
total signal is dominated by the clustering term; the un- 
ambiguous extraction of the single galaxy contribution in 
those cases is difficult if we add another degree of freedom 
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in (7 for each magnitude bin. In the case of m r < 19, the 
best-fit values of a are indeed almost independent of m r 
and a ~ 3'1 as expected from our model assumption. This 
is why we constrained a to be independent of m r in the 
actual fitting procedure. This assumption does not affect 
the best-fit parameters except for E|° with m r > 19. 

3. Implications 

The fitting results presented in the previous section are 
model-independent in a sense that it does not assume any 
a priori relation between the r-band magnitude and the 
FIR emission of galaxies. Therefore the empirical correc- 
tion of the SFD extinction for galaxies can be made from 
Figure 8. We note here that the correction is not limited 
to SDSS galaxies at all, and can be directly applied to 
future galaxy surveys as long as the r-band magnitude of 
galaxies is measured. This is because the dominant clus- 
tering term is contributed by galaxies that are not even 
resolved or identified in the SDSS galaxy catalog. 

Nevertheless it is interesting to consider the un- 
derlying connection between the FIR emission and 
r-band magnitude of galaxies on the basis of the 
present result. Consider a galaxy with intensity profile 

ilOOpm 

extinction should be 



(0)[MJy/sr]. Then its contribution to the r-band 



Y. s g (6;m r 



E(B-V) 



x px I 



(e), 



(10) 



where A r j E(B — V) and p are the conversion factors from 
the color excess E(B — V) to the r-band extinction and 
from 100/im intensity to E(B — V), and are given as 2.751 
and 0.0184, respectively (Schlegel, Finkbeiner & Davis 
1998). Integrating equation (10) over 6 assuming the 
Gaussian profile, we obtain 



xpx / 100Mm [MJy]. (11) 
Finally the 100/im flux, /ioo^m is translated to the 100/im 



27R7 2 £f (m r ) 

nally the 
magnitude: 

mioOMm = -2.51og(/ 100Mm /3.63 x 10- 3 [MJy]), (12) 
and equation (11) is rewritten in terms of m r as 

£|°(m r ) = 36.0x l0- a4mi ™<™ (^-pj . (13) 

Since those magnitudes, m r and mioo/mi, should corre- 
spond to the same galaxy, their difference is equivalent to 
the ratio of their absolute luminosities, Lioo/jm/ L r . Thus 
estimated ratios are plotted in Figure 12. The fact that 
the ratio for a single galaxy is approximately constant in- 
dicates that the statistically averaged ratio of FIR and 
optical luminosities of galaxies are independent of the r- 
band magnitude, which is very reasonable. For compar- 
ison, we plot the ratio for adding the clustering term in 
FIR. In this case the integration of £g(0;m r ) over docs 
not converge because our assumed value of 7(= 0.75) is 
valid only for angular separation less than 1°. Thus we 
evaluate the flux simply by multiplying 2ito~ 2 as in the 
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Fig. 12. Ratio of FIR and r-band luminosities as a function 
of the r-band magnitude of the central galaxy. Crosses indi- 
cate the ratio for single galaxy term, while circles include the 
clustering term as well. 



case of the Gaussian profile. Therefore L 



100/im 



includes 



the contribution of other galaxies, but the total ampli- 
tude is subject to change depending on the more accurate 
profile at larger angular scales. The total ratio follows a 
clear single power-law and we believe that the wiggles of 
the ratio for single galaxies is not real but comes from the 
difficulty in separate the single galaxy contribution from 
the total signal as mentioned at the end of the previous 
section. 

The relation between FIR and optical luminosities of 
galaxies can be directly probed from the sample of 3304 
galaxies overlapped in the SDSS and PSCz (Saunders et 
al. 2000). We select the brightest SDSS galaxy locating 
within 2 arcmin from each PSCz galaxy as its counterpart. 
This procedure enables us to find SDSS optical counter- 
parts for ~ 95% of the PSCz galaxies that are located in 
the SDSS survey region. We apply the K-correction based 
on the "K-corrections calculator" service (Chilingarian, 
Melchior, & Zolotukhin 2010) for r-band, and extrapo- 
late the FIR spectral energy densities using second-order 
polynomials determined from 25, 60 /zm flux for 100 /zm 
(Takeuchi, Yoshikawa & Ishii 2003). Figure 13 is a scatter 
plot of Lioo,an (PSCz) and L r (SDSS) for the PSCz-SDSS 
overlapped sample of galaxies. Figure 13 indicates an ap- 
proximate linear relation between Lioo/mi and L r albeit 
with considerable scatters. The solid lines correspond to 
LiootJ,ni/L r = 20 and 50 for reference. Thus the approxi- 
mate linear relation implied from Figure 12 is largely con- 
sistent with Figure 13 since the IRAS galaxies may be 
preferentially biased toward the FIR brighter ones than 
average. In turn, this confirms that our interpretation 
that £| represents the contribution of a single galaxy. 
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Fig. 13. Relation between L 



lOOfim 



and L r for the 



IRAS/SDSS overlapped galaxies. The solid and dashed lines 
indicate -I/ioO/jm/ir = 20, and 50, respectively. 

4. Stacking analysis of SDSS DR6 quasars 

Y07 did not find any definite anomaly between the 
quasar surface densities and Asfd beyond the statistical 
errors. Nevertheless it is interesting to repeat the stack- 
ing analysis to SDSS quasar sample as well. Indeed as we 
show below, we found that a weaker but similar pattern 
of the enhanced extinction around stacked quasar images. 

For this purpose, we use the SDSS DR6 photomet- 
ric quasar sample (Richards et al. 2009a; Richards et al. 
2009b). The analysis method is basically identical to that 
performed in subsection 2.3 except that we have to use 
the larger magnitude bins (Am r = 1.0) due to the limited 
number of the quasars as well as the weaker signature of 
the effect. The stacked images are plotted in Figure 14. 
As in the case of galaxies, we fit the radial profile to equa- 
tions (1), (6), and (8) assuming that a is independent of 
m r . We find the best-fit value of a = 3.'13, which is similar 
to that for galaxy radial profile. The radial profiles and 
the best-fit curves are plotted in Figure 15. 

Unlike galaxies, the profiles are not completely circular, 
and also the best-fit parameters do not exhibit regular 
behavior as a function of m r . Part of the behavior may 
be due to the contamination of non-quasars objects in 
the photometric quasar sample. Therefore it would be 
better to repeat the analysis for the spectroscopic quasar 
sample, which we plan to do in due course. Nevertheless 
the results indicate a clear signal around the center of all 
the tacked images. If we look at Figure 9 in Y07 carefully, 
a very weak anomaly may be recognized for photometric 
quasars as well. This would be consistent with our current 
finding of the FIR emission around those quasars in the 
SFD map. 

It is interesting to ask if the detected FIR emission 
around quasars originated from (1) quasars themselves, 
(2) their host galaxies, (3) neighbor galaxies due to the 



Fig. 16. Best-fit parameters characterizing the FIR emission 
of quasars against their r-band magnitude. 
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Fig. 17. Ratio of the clustering term and the central quasar 
contribution as a function of the r-band magnitude of the 
central quasar. 

quasar-galaxy and/or quasar host galaxy-galaxy correla- 
tion, and/or (4) some other effects (lensing, for instance). 
Observational studies of quasar V-band luminosity M q 
and that of the host galaxy M g imply a very weak cor- 
relation, at most, with significant scatters. Typically a 
quasar is one or two magnitudes brighter than its host 
galaxy (Hamilton, Casertano, & Turnshek 2008; Letawe, 
Lctawe, & Magain 2010). Given those combined with the 
results for galaxies discussed in §3, the signal may be as- 
cribed to the possibility (3). This deserves further in- 
vestigation, and we are currently working on the analysis 
using the WISE (Wright et al. 2010) and IRAS (Wheelock 
et al. 1994) data to explore the quasar-galaxy correlation 
in FIR. 
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Fig. 15. Radial profiles of stacked quasar images corresponding to Fig. 14. Solid curves indicate the best-fit model of equation (1), 
(6), and (8). 



5. Conclusions 

We have detected the small but systematic contami- 
nation in the SFD Galactic extinction map due to FIR 
emission from galaxies. The amount of the contamination 
of the order of mmag is very small and may be negligible 
for most astronomical purposes. Nevertheless, the con- 
tamination is intrinsically correlated with the large-scale 
structure of the universe traced by galaxies, and therefore 
should be kept in mind in precision cosmological measure- 



ments. Indeed even such a small contamination leads to 
a relatively large anomaly of the galaxy surface number 
density as a function of ^4sfd discovered by Yahata et al. 
(2007). 

Our major result is that a galaxy of r-band magnitude 
m r has an additional contribution to the SFD Galactic 
extinction by an amount of 

AA r (m r ) = 0.087 x io°- 41 ( 18 - m -) [mmag], (14) 

due to the FIR emission from itself (single term), and 
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AA r (m r ) = 0.64 x io°- 17 ( 18 - m -) [mmag], 



(15) 



including the contribution from nearby galaxies (cluster- 
ing term: Figure 18). Note that since the SFD determina- 
tion of conversion factor p has statistical and systematic 
uncertainties of approximately 8%, equation (14) and (15) 
would have the similar level of uncertainties. 

One possible method to correct the SFD map is to sub- 
tract all the contributions of identified galaxies; the cor- 
rected extinction at position 9 is obtained as 



A' r (0) = A r ,s FD (0) - Y^Wi - °' m r^ 



(16) 



where 6j is the position of the j-th galaxy with its r-band 
magnitude of m J r . Note that this correction ignores contri- 
butions of other than the identified galaxies; if the SDSS 
DR7 photometric galaxy catalog is used, only contribu- 
tions from those galaxies listed in the catalog are cor- 
rected. We could use £g 0t (f9j — 0; m 3 r ) in equation (16), 
including the clustering term S|. In that case, however, 
it would over-count the true contribution especially when 
the pixel has many nearby galaxies. In turn it is interest- 
ing to examine the extent of which the anomaly detected 
by Y07 is explained by the SDSS DR7 galaxies alone, by 
applying the correction on the basis of equation (16). 

Most likely improved correction formulae than equa- 
tion (14) and (15) may be obtained by further dividing 
galaxy samples according to their colors, fully exploit- 
ing the multi-band photometries of SDSS DR7. Also 
similar stacking analyses of various current/future sky 
maps in different wavelengths, e.g., WISE, IRAS, AKARI 
(Murakami et al. 2007), WMAP (Bennett et al. 2003), 
Planck (Tauber et al. 2010) among others, are interesting 
to make sure of their validity and/or put constraints on 
the possible contamination. 

We would like to emphasize also that the present 
methodology is useful to derive statistical relations among 
galaxy parameters in different bands and to infer the an- 
gular cross-correlation of galaxies, much beyond the mag- 
nitude limit for individual galaxy surveys. The results for 
those directions will be presented elsewhere. 
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Appendix. IRAS Point Spread Function 

In subsection 2.3, we decompose the stacked radial 
profiles into single and clustering terms, assuming the 
Gaussian PSF. Since the IRAS PSF is known to be very 
complex, we need to check the validity of this assump- 
tion. For the purpose of determining the PSF directly, we 
perform similar stacking analysis with SDSS stars. 

More specifically, we select spectroscopic stars brighter 
than m r = 17.0 (12823 stars in total) from the SDSS DR7 
catalog. We first stack the SFD map centered on those 
SDSS stars as we did for SDSS galaxies and quasars, but 
we find no significant signature. This is mainly because 
the bright point sources in the IRAS catalogue are al- 
ready removed from the SFD map. Therefore we go back 
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Fig. 19. Radial profiles of stacked star images. The dashed 
curve indicates the best-fit Gaussian profile. 

to the original ISSA 100/im diffuse map, and perform the 
stacking analysis. The resulting stacked average radial 
profile of the SDSS stars is shown in Figure 19. The data 
points are well approximated by a single Gaussian. The 
best-fit Gaussian shown in the dashed curve has a width 
of a = 2'A2, which is slightly smaller than that found 
in subsection 2.3. This is understood because since the 
SFD map is constructed by further smoothing the origi- 
nal ISSA map. Thus we conclude that our assumption of 
the Gaussian PSF is valid. 
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